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Conformal Weyl gravity and perihelion precession 
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We investigate the perihelion shift of planetary motion in conformal Weyl gravity using the metric of 
the static, spherically symmetric solution discovered by Mannheim and Kazanas [Astrophys. J. 342, 635 
(1989)]. To this end we employ a procedure similar to that used by Weinberg for the Schwarzschild 
solution, which has also been used recently to study the solar system effects of the cosmological constant 
A. We show that besides the general relativistic terms obtained earlier from the Schwarzschild-de Sitter 
solution, the expression for the perihelion shift includes a negative contribution which arises from the 
linear term yr in the metric. Using data for perihelion shift observations, we obtain constraints on the 
value of the constant y similar to that obtained earlier using galactic rotational curves. 
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I. INTRODUCTION 

One of the possible alternatives to the standard second- 
order Einstein theory of gravity which has been proposed 
during the last two decades is conformal Weyl gravity 
[1—3]. Instead of fixing the gravitational action by demand- 
ing that the theory be no higher than second order as in the 
case of the Einstein-Hilbert action, Weyl gravity employs 
the principle of local conformal invariance of spacetime as 
the supplementary condition that fixes the gravitational 
action. This means that the theory is invariant under local 
conformal stretching of spacetime geometry of the form 

gp.v(x) -» Cl 2 (x)g /L1 ,(x), (1) 

where fl(x) is a smooth, strictly positive function. This 
restrictive conformal invariance leads to a fourth-order 
theory with a unique, conformally invariant action 

I w =-«f d 4 x(-g) l/2 C^ VK C^ VK 

= -2a f d 4 x(-g) 1 / 2 [R flK R^ K - (R^) 2 / 3] 

+ a total derivative, (2) 

where C A(LtJ/K is the conformal Weyl tensor and a is a 
purely dimensionless coefficient. This action gives rise 
[1] to the gravitational field equations given by 

yFUg^agvp 4^ = -2 aW^ v = - ' T IAl „ (3) 

°gct(3 Z 
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where T /J t , is the stress-energy tensor and 

= 2 C\/. ap + C%/R a p, (4) 

is the Bach tensor. When R jJV is zero, the tensor W ^ v , 
which consists of R^ v and its derivatives, vanishes, too, so 
that any vacuum solution of Einstein’s field equations is 
also a vacuum solution of Weyl gravity. However, the 
converse is not necessarily true. Despite the highly non- 
linear character of the field equations, a number of exact 
solutions [4-7] of conformal Weyl gravity has been found 
in the case of spherical and axial symmetry. Recently, 
cylindrically symmetric solutions [8,9] in Weyl gravity 
have also been studied. Moreover, by studying the interior 
structure of a static, spherically symmetric gravitational 
source, it was shown [7] that the field inside the source is 
described exactly by a fourth-order Poisson equation, 
which under the proper conditions still admits a 
Newtonian potential \/r term. Therefore, although the 
second-order Poisson equation in general relativity is suf- 
ficient to generate a Newtonian potential, it is not by any 
means a necessary requirement, so that Newton’s law of 
gravity remains valid in the fourth-order Weyl gravity. 

The exact static and spherically symmetric vacuum 
solution for conformal gravity is given, up to a conformal 
factor, by the metric [1] 

d r 2 

ds 2 = -B{r)dt 2 + — + r 2 (dd 2 + sin 2 0d(f 2 ), (5) 

where 

B(r) = 1 ~ 13(2 - 3/3y> - 3 (3y + yr - kr 2 , (6) 
r 
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and (3, y, and k are integration constants. This solution 
includes as special cases the Schwarzschild solution 
(y = k = 0) and the Schwarzschild-de Sitter (y = 0) so- 
lution, the latter requiring the presence of a cosmological 
constant in Einstein gravity. The constant y has dimensions 
of acceleration and so the solution provides a character- 
istic, constant acceleration without having introduced one 
at the Lagrangian (such as in the relativistic implementa- 
tion of MOND with TeVeS [10]). It should be noted that a 
metric with a similar linear term was also introduced 
recently [11] in a model for gravity at large distances 
based, however, on very different considerations. 

Although the magnitude and the nature of the integration 
constant y in (6) remain uncertain, it has been associated 
[1] with the inverse Hubble length, i.e., y — 1 /R H . In this 
case the effects of the acceleration are comparable to those 
due to the Newtonian potential term 2(3/ r = rj r ( r s is the 
Schwarzschild radius) on length scales given by 


rjr 2 - y - \/R„ or r - {r s R H ) x '\ (7) 


As noted in Ref. [1], for a galaxy of mass M — 10 11 M G with 
r s ^ 10 16 cm and R H ^ 10 28 cm, this scale is r ~ 10 22 cm, 
i.e., roughly the size of the galaxy, a fact that prompted 
Mannheim and Kazanas to produce fits to the galactic 
rotation curves using the metric of Eq. (6) above. One 
should point out that Eq. (7) describes not a particular length 
scale but a continuum of sizes at which the contribution 
from the linear term becomes significant. Objects along this 
sequence encompass not only galaxies but also, at larger 
scales, galaxy clusters and, at lower scales, globular clus- 
ters, only recently found to require the presence of dark 
matter in order to account for the observed dynamics 
[12,13]. As we know, to account for the dynamics of objects 
along this sequence within the standard gravitational theory, 
one would need to invoke the presence of dark matter. 

The classical tests in Einstein’s gravity for the metric (6) 
with vanishing y have been well documented since the 
discovery of the theory and, therefore, it would be natural 
to study any effects of the linear term in the metric on such 
tests. The issue of the propagation of null geodesics, par- 
ticularly the computation of the bending of light by a 
spherically symmetric object using the metric in (6), has 
been studied in detail. The first results obtained in 
Refs. [14-17] showed that besides the positive Einstein 
deflection of 4(3 /b, the expression for the deflection of 
light in Weyl gravity contained an extra term — yb, where 
b is the impact parameter. This led to the paradoxical 
situation where the bending angle in lensing increased 
with the light ray’s impact parameter with respect to the 
lens. This problem was later solved in Ref. [18], where it 
was shown that when the curvature of the background, 
asymptotically nonflat geometry of (6) is taken into 
account, the total bending angle is given by 


A iff 


4/3 _ 2(3 2 y _ kb 3 
~b~ b ~ 2/3 ’ 


(B) 


such that the contribution from the linear term in the metric 
is inversely proportional to the impact parameter. 
Moreover, its ratio to that of the standard 1/r component 
is of order (3y, which, given the associations and magni- 
tudes of these constants for a galaxy, we get (3y — 10“ 12 , 
i.e., insignificant for all practical purposes. 

In this paper we study timelike geodesics in Weyl grav- 
ity, particularly the effect of the linear term in the metric on 
the perihelion shift. Then using the available data for 
planetary perihelion shifts, we get constraints on the mag- 
nitude of y similar to that obtained earlier [1] from the 
fitting of galactic rotational curves. Earlier studies [19-22] 
of the geodesic structure in the Schwarzschild-de Sitter 
solution showed that the cosmological constant A in- 
creases the perihelion shift in the Schwarzschild solution 
by 77 A 7 ( 1 — e 2 ) 3 /m, where a, e are the length of the 
semimajor axis and eccentricity, respectively, and m = 
GM/c 2 , with M being the mass of the source. Some 
authors [23,24] have used this correction to the perihelion 
shift of Mercury to obtain upper bounds for the cosmologi- 
cal constant, while others [25,26] showed that the effect of 
the cosmological constant is only significant at large radii 
and since it is unmeasurably small for Mercury’s orbit, it 
cannot be used to limit the cosmological constant. 

In Sec. II we start with the differential equation repre- 
senting timelike geodesics and use a method similar to that 
used by Weinberg for the Schwarzschild solution to obtain 
an expression for the perihelion shift of a test particle in the 
exterior geometry described by the metric in (6). Then in 
Sec. Ill, the results are summarized and discussed. 


II. TIMELIKE GEODESICS, 
PERIHELION PRECESSION 

The timelike orbits for the metric 

ds 2 = - B{r)dt 2 + A(r)dr 2 + 7(dd 2 + sin 2 ddcf> 2 ), (9) 


are given by [27] 

A{r)( dr \ 2 1 _ 1 _ _ E 

~d r \d(f)) 7 J 2 B{r ) ~~ T 2 ’ 


( 10 ) 


where J and E are constants of the motion. Since (3y <C 1, 
for simplicity we take 

2(3 

B(r) = A~ l (r) = 1 - — + yr - kr 2 . (11) 


Then the angular distance between the perihelion r_ and 
aphelion r + is given by 




1 l-t/2 dr 

7\ 7' 


( 12 ) 
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Using the fact that dr/d<f> vanishes at r_ and r+, one can derive the following values for the constants of the motion 

A, A 


E = 


B(r+) BOA 


(13) 


and 


j2 _ B(r + ) B(r_) 

J 1 _ 1 

^ 7T 


such that the angular distance in (12) becomes 


<t>(r+) ~ 4>{r~) = J + A(r) 1/2 ^ 


r±(B~ l (r) ~ B~ l {rJ)) - r 2 + (B~ x (r) - £'V+)) 1 1 ~ l / 2 dr 


r 2 + rl (B“ 1 (r + )-5“ 1 (r_)) 


ji ■ 


(14) 


(15) 


The perihelion precession per orbit is given by 

A(f> = 2\<f>(r+) - </>(r_)| - 277. (16) 

The expression in the square brackets in (15) vanishes at r + and r_ and, hence, for slightly eccentric orbits we can write 
\A{B-\r) - - r 2 + (B~\r ) - B~\r+)) 1 


L r 2 + r 2 L{B-\r+) - B~\r_)) 

where C is a constant. Then letting u = \/r and differ- 
entiating twice with respect to u gives 




(17) 


<f>(r+) ~ 4>(r-) = - f 

J U- 


A(u,y/ 2 du 

[C(w_ — u){u — M + )] 1 / 2 


C * 1 - 


(u + — m_)(m + + u.-)A"(u ) 
2 (A(u+) - A(u-)) 


Using the substitution 


u=L 


where L = 2/(u+ + m_) = a( I — e 2 ) is the semilatus 
rectum of the elliptic orbit. Again, for slightly elliptic simplifies the integral to 
orbits we can write 


1 . 1 

u = -(m + + u ~ 1 -(«+ — M_)stni/^, 


so that 


. . . ,/u+ + U-\ 

A(u+) - A(u_) ~ (; u+ - uJ)A'y j, 

uA"(u) 


C « 1 - 


A'(u) 


(19) 


( 20 ) 


<f>(r+) ~ 4>(r-) = — J ' ^A{ipy/ 2 d<p. 


This leads to 


<f>(r+) ~ ) ~ 77^1 + 


— 7t/2 

3/3 


(24) 


(25) 


(26) 


17 


a( 1 - e 1 ) 


2 +^ka 2 (l-e 2 ) 2 


Now 

A(u) = B~ l (u) 

= 1 + 2 (3u + 4 0 2 u 2 + — - - + X ~ 4 (3y + 

U ll u 2 


Hence, substituting this in (20) we get 

2/3 + 3 yu~ 2 — 12 k(3u~ 2 — 8 ku~ 2 


3k 

w 


+ ^lfl 3 (l - e 2 ) 3 - — fl 2 (l - e 2 ) 2 


— 2ya(l — e 


L) > 


2/3 


(27) 


( 21 ) 


C 


2/3 + 8 (3"u — 4k(3u 2 + yu 2 — 2 ku 


i _ 3 


u=L 1 

( 22 ) 


so that the perihelion shift per orbit is given by 

A </> ~ + ^ka 2 ( 1 - e 2 ) 3 - %« 2 ( 1 - e 2 ) 2 . 

a( 1 - e ) (3 (3 

(28) 


III. DISCUSSION AND CONCLUSION 


or 


The expression for the perihelion precession in confor- 
mal Weyl gravity obtained in (28) constitutes the main 
result of this paper. Besides the conventional term 
(23) 6n(3/a(l — e 2 ) representing the precession in the 

Schwarzschild geometry, the expression includes two other 
The expression for the angular distance in (15) can be terms that originate due to the large-scale structure of the 
written as embedding spacetime. The /(-dependent term is clearly 


2y y 4k 3k 

C ~ l — 4 /3m — — + — =■ + — — — r 
m /3u- u 2 /3u J 


u=L 
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associated with the kr 2 de Sitter term of the metric of Weyl 
gravity and is the same contribution arising from the 
cosmological constant obtained earlier in Refs. [19-22] 
for the Schwarzschild-de Sitter solution. This contribution 
was expected considering the fact that the metric in (6) 
reduces to the Schwarzschild-de Sitter solution when y 
approaches zero. The additional term Trycr(\ — e 2 ) 2 / (3 in 
(28) originates from the linear yr term in the metric, and 
due to its negative sign, it reduces the amount of precession 
per orbit. This has also been observed before in Ref. [18] 
[see Eq. (8) in this paper], where it was shown that both the 
linear and de Sitter terms in the metric diminish the light- 
bending angle. In this case, however, the effects on the 
precession from these two terms are opposite. 

The nature and magnitude of the constant y are 
still unknown. When (3 = 0 or when r is sufficiently large 
so that all (3— dependent terms in (6) can be ignored, 
the metric can be rewritten under the coordinate trans- 
formation 


of a spherically symmetric matter distribution in a cosmo- 
logical background. Assuming a spatially flat matter- 
dominated universe with k = 0, one can use (28) to get 
an upper bound for the magnitude of y from the difference 
between the observed and the general relativistic values of 
the precession of perihelia. So, for example, in the case of 
Mercury [28], 

A</> obs — A (pgr = —0.0036 ± 0.0050 arcsec/century, 

(31) 


where A </> ohs refers to the observed value of the perihelion 
precession once corrected for the general precession of the 
equinoxes and for the perturbations due to other planets. 
The contribution to the precession of the perihelion from 
the linear term in the metric can be written in the form 


A 4>y 


1 77 y / 

r s 


(32) 


2(1 + yr 
j R(t)dt, 


4 r 

kr 2 ) 1 / 2 + 2 + yr’ 


and 


in the form 


(29) 


ds 2 = 


[1 - p 2 (y 2 /16 + k/ 4)] 2 
* 2 (t)[(1 - yp/4) 2 + kp 2 / 4] 2 

r , 7 r 2 (t) 

X -dr 2 + 


[ “ ' ' [1 ~ p 2 (y 2 /16 + k/ 4)] 2 

X (dp 2 + p 2 (d0 2 + sin 2 0d<f) 2 )) j. 


(30) 


The metric is, therefore, asymptotically conformal to a 
FLRW metric with arbitrary scale factor R(t) and spatial 
curvature k = — k — y 2 /4, i.e., it describes a spherically 
symmetric object embedded in a conformally flat back- 
ground space. The fact that the curvature of this background 
space depends on y and k points toward a cosmological 
origin of y. This led to its association with the inverse 
Hubble length in Ref. [1], so that y ~ 10 -28 cm -1 , and 
this explained quite successfully the flat rotation curves in 
galaxies and galaxy clusters. 

On the other hand, there is nothing in the theory that 
forbids y from being also system dependent (i.e., like the 
parameter (3 , which is associated with the mass of the 
gravitating body), and in this case one can suggest that 
the linear term yr in the metric provides the necessary 
changes in the spacetime geometry to allow the embedding 


where r s = 2(3 = For Mercury r s L 1 = 5.3 X 10 8 
so that 

A <py = 1.2 X 10 8 yL rad/orbit. (33) 

Therefore from (31) we get y < 1.5 X 10 -31 cm -1 , which 
is about three orders of magnitude smaller than the large- 
scale estimate obtained from the fitting of galactic rota- 
tional curves. A similar calculation for other planets such 
as the Earth leads to tighter upper bounds for y. However, 
when comparing bounds on y derived from planetary data 
with those obtained from the large-scale geometry of the 
universe, one should keep in mind that the calculation of 
the perihelion precession assumes that planets are treated 
as test particles devoid of self-gravity, and so this weakens 
the validity of such bounds. In fact, in the case of the 
asteroid Icarus with a diameter of just 1 km, (so that 
the effect of self-gravity in this case is insignificant), the 
difference between the observed and general relativistic 
perihelion precession is given by [27] 

A0 obs — A^> gr = —0.5 ± 0.8 arcsec/century, (34) 

which gives, using (32), the value y < 1.3 X 10“ 28 cm -1 . 
This fits nicely with the cosmological estimate of 
10 -28 cm -1 obtained earlier from galactic rotational 
curves. 
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